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Показано, что в зависимости от создаваемых условий тепломассообмена в пористом слое древесных отходов реализуются низ- 
ко- и высокотемпературный режимы горения, и режим тления, при котором основная масса газообразных компонентов, обра- 
зующихся при пиролизе, выбрасывается в атмосферу, не сгорая. С помощью математической модели, рассмотренной в работе, 
можно прогнозировать режимы сжигания пористых, коксующихся отходов и определять состав образующихся при этом газооб- 
разных продуктов. 


Самым распространенным способом утилиза- 
ции строительных отходов является их сжигание в 
волне горения [1,2]. Проблема утилизации отходов 
с оптимальными экологическими последствиями 
является актуальной и сводится к выбору тепло- 
массообменных режимов воздействия на топливо, 
чтобы количество вредных газов в продуктах сгора- 
ния было минимальным. 

В данной работе рассматриваются процессы 
утилизации отходов деревообработки в вертикаль- 
ной цилиндрической печи заданных размеров. 

Физико-математическая модель. Предположим, 
что твердое пористое топливо (отходы деревообра- 
ботки) помещено в расположенную вертикально 
цилиндрическую печь. Обозначим высоту слоя то- 
плива через И. После загрузки осуществляется про- 
дувка печи высокоэнтальпийным газовым потоком. 

Будем считать, что процессы сушки и пиролиза 
твердого горючего описываются уравнениями вида: 
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Здесь ѵ,Д/| - масса древесных отходов, ѵ 2 М 2 - 
воды, ѵ\М\ - водяного пара; ѵ\М 2 - кокса, ѵ\М 4 - 
золы, ѵ\ а М І0 - «-компоненты газовой фазы, обра- 
зующейся при пиролизе древесины. 


Кроме данных реакций учитывались гетероген- 
ная реакция горения кокса и гомогенные реакции 
окисления оксида углерода, водорода и метана: 

С + О, = СО, + <7 С , 2С0+0, = 2СО, + <7 С0 , 
2Н,+0, =2Н,0 + <7 н , 

СН 4 + 20 2 — СО, + 2Н,0 т- <у сн . 

В приведенных выше формулах </,, ^ Р , д с , ^ со , <? Н2 , 
ц сн - тепловые эффекты реакций испарения воды, 
пиролиза древесины, окисления кокса, оксида 
углерода, водорода и метана соответственно. 

В силу цилиндрической симметрии физической 
задачи математическая модель, описывающая дан- 
ный процесс, была записана в цилиндрической си- 
стеме координат. Теплообмен на боковой поверх- 
ности на входном и выходном сечении рассчитыва- 
лись с использованием граничных условий третье- 
го рода. Предполагая, что изменение искомых 
функций по радиусу значительно меньше, чем 
вдоль оси печи уравнения осреднялись по радиусу, 
в результате была получена следующая система 
уравнений: 

• сохранения массы исходного конденсирован- 
ного вещества 

Р '^1~ = ~ Кр ' К г =к г р &' і ех р(“ІУ’ (1) 
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• сохранения массы воды в связанном состоянии 

Рі ууу = ~ к ь > = к ь р 2 (р 2 ехр -^-1 , (2) 

• сохранения массы кокса 

д ю, 

Рз ^ ^ — сс с К р — $К С , 

К = ^ Р5<Р5 с А<Рз ехр -Ц- | , (3) 


1 = 0- г ' = 1 ’ 4 > с а =с аЯ , а =1,6, 
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( 10 ) 


2 = 0:7’ = Т Е ,р = р Е ,Х 8 ^ = а 1 (Т,-Т Е ),с а =с аЕ , 

, д с а дТ д Т :і 
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В уравнениях (6-8) введены краткие обозначе- 
ния скоростей реакций 


• сохранения массы золы 
Ра 
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неразрывности газовой фазы 
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сохранения массы а-компоненты газовой фазы 
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• сохранения энергии газовой и конденсирован- 
ной фазы 
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движения газа в порах (в форме квадратичного 
закона фильтрации) и уравнение состояния 
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Система уравнений (1-9) решалась при следую- 
щих краевых условиях: 
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Здесь ѵ - скорость фильтрации; а в , а с , а 5 - до- 
ля газообразных продуктов, кокса и золы при пиро- 
лизе; М 5 - молекулярная масса газовой фазы; ин- 
декс 5 приписывается величинам, характеризую- 
щим газовую, а с- конденсированную фазу; к - ко- 
эффициент проницаемости пористого слоя; Т 5 , 
Т - температура конденсированной фазы и газа в 
порах; Е Р , Е ь , Е с , Е со , Ещ, Е С щ, к Р , к ь , к с , к со , кщ, к С щ 
- энергии активации и предэкспоненциальные 
множители реакций пиролиза древесины, испаре- 
ния воды гетерогенной и гомогенных реакций оки- 
сления С, СО, Н 2 и СН 4 соответственно; К - уни- 
версальная газовая постоянная; і) а , Я 5й , с а , М іа - ко- 
эффициент молекулярной диффузии, теплопро- 
водности, массовая концентрация и молекулярная 
масса а-компоненты, параметр с индексом а=1 со- 
ответствует 0 2 , а=2 - СО, а=3 - СО,, а= 4 - Н 2 , 
а= 5 - СН 4 , а=6 - Н,0, а=7 — 1Ч 2 ; р — давление га- 
зообразных продуктов в порах; Р=р/р н - безразмер- 
ное давление; р н - характерное давление; А х =а Р ? - 
объемный коэффициент теплоотдачи; - удельная 
поверхность пор; ^ К =е 8 аТ^ - радиационный поток 
от внешней стенки печи в окружающую среду; 
е 5 — степень черноты; 1— время; у = 7 / 1 *— безразмер- 
ная пространственная переменная, і - размерная 
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пространственная переменная, направленная по 
оси печи; а ь а 2 - коэффициенты теплоотдачи, ха- 
рактеризующие конвективный теплообмен во вход- 
ном (выходном) сечении и на внешней стенке печи; 

/, = у](\КТ, 2 ехр(Е С0 1 КТ,)І(р х к с ^ со Е с0 )) - мас- 
штаб длины; Т, - характерная температура; 
<р 0 =(/Т 2 -г 0 2 )/г 0 2 , г х - наружный, г 0 - внутренний ради- 
ус печи; р 0 , с р „, - плотность, теплоемкость и тепло- 

проводность стенки печи; р„ с рі , А„ <р, - плотность, 
теплоемкость, коэффициент теплопроводности и 
объемная доля древесины (г-1), воды (/=2), кокса 
(г-3), золы (г-4) и газовой фазы (/=5); Я 5 , Я й , Я^- ко- 
эффициенты молекулярной и лучистой теплопро- 
водности газовой фазы и коэффициент теплопро- 
водности конденсированной фазы; а - постоянная 
Стефана-Больцмана; р - коэффициент динамиче- 
ской вязкости; р - константа в квадратичном законе 
фильтрации; # - ускорение свободного падения; 
и=ѵ/ѵ*~ безразмерная скорость, масштаб скорости 
ѵ*=0, 797- КИм/с; масштаб длины /*=0,95-1(Им. 

Краевая задача (1-11) решалась численно. Тех- 
нология решения задачи вида (1-1 1) с использова- 
нием итерационно-интерполяционного метода [3] 


приведена в работах [4, 5]. Кинетические параме- 
тры процессов сушки, пиролиза и горения, а также 
теплофизические коэффициенты газообразных 
компонентов и конденсированных фракций были 
взяты из работ [6-13]. 

Результаты расчетов и их анализ. При численном 
решении задачи были найдены низко- и высокотем- 
пературный режимы горения, и режим тления. Эти 
режимы зависят от условий тепломассообмена по- 
ристого слоя с продуваемым через него окислите- 
лем. Низкотемпературный режим реализуется при 
слабом вдуве не очень горячего окислителя. Профи- 
ли температуры и объемной доли влаги при низко- 
температурном горении для шести моментов време- 
ни, представлены на рис. 1 . Данный режим горения 
был получен при 7^=800 К и АР=Р Е -Р Н =0,005 
(Др=500 Па). Здесь кривые 1-3 рис. 1, а, соответ- 
ствуют процессу сушки, а кривые 4-6 процессу об- 
разования и распространения волны горения. Из 
данного рисунка следует, что при рассматриваемых 
условиях зажигание происходит на тыльной стороне 
слоя, а затем волна горения распространяется про- 
тив потока. На рис. 1, 6, приведены профили объе- 
мных долей воды для тех же моментов времени. 




Рис. 1. Профили температуры: а) газовой фазы и б) объемной доли влаги в моменты времени I, с: 1) 4,8 ; 2) 56,9; 3) 145,1; 
4) 311,7; 5) 390,3; 6) 482,1 




Рис. 2. Профили: а) объемной доли отходов деревообработки и б) скорости фильтрации в моменты времени I, соответствую- 
щие рис. 1 
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Анализируя данный рисунок, видим, что к мо- 
менту времени соответствующему кривой 4 
рис. 1, а, объемная доля влаги равна нулю во всем 
объеме, после чего происходит зажигание. 

Профили объемной доли исходного конденси- 
рованного вещества приведены на рис. 2, а, кото- 
рый подтверждает описанный выше механизм за- 
жигания и горения отходов деревообработки. Здесь 
кривые 1-5 соответствуют кривым 1, 3-6, рис. 1, а. 
Видим, что топливо зажигается с тыльной стороны 
(т=/г), и горение распространяется в направлении, 
противоположном направлению вдува горячего 
окислителя. На рис. 2, б, изображены профили 
скорости фильтрации для моментов времени, соот- 
ветствующих кривым 1-3 и 6, рис. 1, а. В зоне ис- 
парения скорость фильтрации имеет локальный 
максимум (кривые 1-3), а затем после высушива- 
ния слоя профиль скорости становится монотонно 
возрастающим (кривая 4, рис. 2, б). 

Из анализа кривых 1-4, рис. 3, а, которые соот- 
ветствуют кривым 1-4, рис. 1, а, можно заключить, 
что к моменту зажигания в порах топлива пары во- 
ды отсутствуют. 


На рис. 3, б представлены профили массовой кон- 
центрации в моменты времени, соответствующие 
кривым 1, 3 и 6, рис. 1, а, т. е. кривые 1 и 2 соответ- 
ствуют процессу сушки, а кривая 3 - моменту рас- 
пространения волны низкотемпературного горения. 
При перепаде давления АР=0,005 и температуре 
Г г =1000 К реализуется также низкотемпературный 
режим горения. Но в отличие от предыдущего низко- 
температурного режима волна горения распростра- 
няется в направлении движения потока окислителя. 

Высокотемпературный режим горения реализо- 
вывается при более высокой скорости продувки, 
например, взяв ДР=0,05, а температуру вдуваемого 
окислителя Т Е - 1000 К. На рис. 4, а, приведено 
пространственное распределение температуры при 
данном режиме горения для четырех моментов вре- 
мени: прогрева (кривая 1), образования (кривая 2) 
и распространения (кривые 3, 4) волны горения. 
Пространственные распределения перепада давле- 
ния для данного режима горения изображены на 
рис. 4, б. Здесь кривая 1 соответствует стадии про- 
грева, 2 - моменту образования, 3 - началу распро- 
странения волны горения. 




Рис. 3. Пространственные распределения: а) массовых концентраций паров воды и б) диоксида углерода в моменты времени, 
соответствующие рис. 1, а 

Т, К а) Д Р 




Рис. 4. Профили: а) температуры и б) давления для моментов времени I, с: 1) 13,9; 2) 54,4; 3) 81,7; 4) 139,2 
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Профили массовой концентраций СО для вы- 
сокотемпературного режима приведены на рис. 5, 
а. Здесь кривые 1, 2 соответствуют стадиям прогре- 
ва и образования, а кривые 3, 4 моментам распро- 
странения волны горения. В силу того, что профи- 
ли массовых концентраций остальных компонен- 
тов имеют качественно такой же вид, то они не 
приведены. Изменения объемной доли кокса, об- 
разующегося при прогреве, воспламенении и горе- 
нии, представлены на рис. 5, б. 

Представляет интерес режим тления, при кото- 
ром происходит сушка и пиролиз твердого топли- 
ва, а затем окисление коксового остатка вплоть до 
золы, при этом температура в сжигаемом слое не 
превышает температуры вдуваемого окислителя. 
Данный режим реализуется, например, если 
Р Е -Р Н = 0,1 (Л/;=10 4 Па) и Г г =800 К, т. е. при интен- 
сивном вдуве не очень горячего окислителя. На 
рис. 6, а, приведены профили температуры в раз- 
ные моменты времени. Здесь кривые 1-3 соответ- 
ствуют стадии сушки и пиролиза, а кривая 5 соот- 
ветствует моменту времени, начиная с которого 
твердое топливо полностью высушено. На рис. 6, б, 


показано, что при тлеющем режиме сжигания не 
только идет процесс сушки и пиролиза древесных 
отходов с образованием кокса, но происходит и 
окисление кокса. Здесь кривые 1-3 соответствуют 
стадии прогрева и начала пиролиза. Затем (см. кри- 
вые 4 и 5, рис. 6, б) в окрестности входного сечения 
печи начинается реакция низкотемпературного 
окисления древесного угля, и объемная доля кокса 
начинает уменьшаться. 

На рис. 7, а, изображены профили скорости 
фильтрации, из которых следует, что профиль скоро- 
сти, пока не испарилась влага, является нелиней- 
ным. Здесь кривые 1-4 соответствуют процессу суш- 
ки (кривым 1-4, рис. 6, а), а кривая 5 соответствует 
кривой 6, рис. 6, а. И, наконец, на рис. 7, б, приведе- 
ны профили давления вдоль оси реактора. Кривые 1 
и 2 данного рисунка соответствуют начальной ста- 
дии процесса, т. е. кривым 1 и 2, рис. 6, а, кривая же 
3 соответствует кривой 5, рис. 6, а. Следует отметить, 
что пока идет процесс сушки и интенсивное выделе- 
ние водяного пара, профиль давления также нели- 
нейный и становится линейным только после просу- 
шивания пористого слоя топлива. 


5) 




Рис. 5. Профили: а) массовой концентрации оксида углерода и б) объемной доли кокса в моменты времени I, с: 1) 13,9 ; 
2) 54,4; 3) 81,7; 4) 139,2 




Рис. 6. Профили: а) температуры и 6) объемной доли кокса в моменты времени {, с: 1) 1,2; 2) 5, 1; 3) 11,4; 4) 21,9; 5) 79,3; 6) 747 
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Рис. 7. Пространственное распределение: а) скорости фильтрации и 6) давления в порах топлива для режима тления в момен- 
ты времени, соответствующие рис. 6, а 


Следует отметить, что достоверность численного 
алгоритма используемого в данной работе исследо- 
валась ранее в [5], где в качестве пористого коксую- 
щегося топлива брался торф. Численные расчеты 
работы [5] хорошо согласовывались по температуре 
и скорости горения торфа с экспериментальными 
данными работы [14]. Остальные характеристики в 
экспериментальной работе не определялись. 

Выводы 

• Показана возможность реализации различных 
режимов воспламенения и горения отходов де- 
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